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ABSTRACT 


The use of superconductors in the field Windings of large 
Synchronous rotating electrical machines virtually eliminates the 
need for magnetic iron core material and hence reduces the size and 
weight requirements for a given power level. This possible reduction 
in size and weight can be exploited in the propulsion of hydrofoil 
craft. Present hydrofoil propulsion schemes all have serious salient 
disadvantages that have stimulated the search for a better propulsion 
scheme for hydrofoil craft. This thesis proposes the use of a novel 
dual armature superconducting induction motor. The motor would be 
encased in a water tight pod and placed on the end of a foil Strue 
and directly drive a Super=cavitating proreller. A mathematical 
model is developed and the size and parameters of a 21000 horsepower 
motor are determined. The resulting model is statically and 
dynamically analyzed. 
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I INTRODUCTION 
In conventional displacement type ship hull forms, the power availe- 
able within the hull can be readily converted to propulsive thrust by 
conventional direct drive propellers. However, in a hydrofoil type ship 
hull form, the mechanism of converting the available power within the hull 
to propulsive thrust is no longer a simple matter. When the hydrofoil 
vehicle ‘flies’ (becomes foil borne), the hull is lifted out of the water 
and conventional types of propulsion are no longer feasible. A Simple 
direct drive straight rotating shaft with a propeller on the end will 
not suffice. 
An assortment of alternative hydrofoil propulsion schemes have 
been proposed or developed: 
1. The water jet system scoops up water at an inlet in a forward 
foil strut and ejects this water at high velocity out the stern, 
thus providing the propulsive thrust. This arrangement requires 
large amount of internal hull volume for water ducts and pump 
machinery and places large amounts of water weight inside the hull. 
Additionally, this method has a low efficiency at low velocities. 
Ze Another scheme utilizes a mechanical gear train containing 
right angle gears to drive propellers mounted on the foil struts. 
This scheme suffers from a very high initial cost and a low mean 
time between failure, hence unreliability and high maintenance costs. 
36 Air propellers, turbo-props, and other aerodynamic propulsors 
all have los efficiencies in the speed ranges anticipated for 


hydrofoils in the near future. 





4, Electric motors can be placed in pods on the end of the foil 
struts and drive propellers directly. This method eliminates the 
Objections of the preceding schemes, but introduces a new 
objection of its own. The pod of the size required for a conventional 
motor would produce extremely high levels of hydrodynamic drag in 
the speed ranges of the hydrofoil vehicles. 
While all of the above methods of propulsion are feasible and can be 
utilized, there is much needed improvements before hydrofoil vehicles 
can compete economically with the displacement type vessels. Much 
research and development is being done to improve on each type of propule 
Sion system. In particular, the motor in the submerged pod system has 


been vigorously pursued with emphasis on reducing the required size of the 
(4559657) 


The first attempt at applying super conducting materials to motors 


(7) 


motor through the use of super conduction technology. 


utilized Faraday's homopolar concept. It has the characteristics of 
low DC voltages with high current hensitice: ape The high DC current 
densities of the homopolar devices complicates the power transmission and 
switching problems and requires some sort of liquid brushes in the motor. 
several different configurations using the homopolar concept are 
discussed in reference (4). 

The recent trend in super conducting machines has been toward the 
use of AC synchronous devices. While AC machines lack the high starting 
torques which are characteristic of DC machines, they have a superior 
reliability in a marine environment and the power transmission and 


Switching problems are much less than for DC machines. The AC super 


conducting machines most considered are of the type with a stationary three 





V4 
phase armature (stator) and a revolving DC field (rotor) consisting of 


super=-conducting winding in a Dewar bottle. (49597) 

A recent developnent in AC supereconducting machines is a family of 
machines which utilizes two armatures and a supereconducting field 
winding. (2) The field winding and one armature, which may or may not be 
wound, rotate, and the other armature is fixed (stator). One motor config- 
uration utilizes a superconducting field which is free to rotate and is 
mechanically unconstrained except for its own inertia, it responds only 
to the applied electrical torques and its inertia. Both armature windings 
are wound in a three phase manner, one armature is fixed (stator) and 
receives the input electrical power, and the other, positioned between the 
stationary armature and the rotating field, develops the mechanical output 
torque in the same manner as @ wound rotor induction motor does. This 
latter configuration is the topic of this thesis. It will be mathematically 
developed and applied to the propulsion of a hydrofoil vehicle. Its static 
and dynamic characteristics will be presented. The resulting size of the 
motor will demonstrate that it is compact and will significantly reduce the 


hydrodynamic drag of the pod to an acceptable level. 
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II PHYSICAL DESCRIFTION OF THE DUAL ARMATURE MOTOR 

The motor to be studied by this thesis utilizes a four pole 
rotating field that is mechanically unconstrained, except for its own 
moment of inertia. Tnis field will be of the super conducting type and 
will be encased in a Vewar bottle and cooled by liquid helium. It will 
be excited by a DC current sourcee In general, the design of the field 
will be similar to the field winding developed at mrt, (13) This field 
Will be capable of producing conventionally unachievable magnetic flux 
densities due to the large current densities that are made possible by 
the super conducting material. The superconducting winding will hence 
forth be referred to as the field. 

surrounding the field and attached to it is a damper shell. This 
shell will nrovide the start up toraves and damp out dynamic oscillations. 
The damper shell and the field together will hence forth be referred to as 
the rotor. 

surrounding the rotor will be a three phase, four pole winding. 
This winding will contain no magnetic core material and will contain 
only enough steel or other material to provide mechanical support and 
rigidity. it will be free to rotate relative to the rotor and to the 
motor frame. It will be attached to a shaft from which mechanical power 
may be extracted and supplied to the load, in this case a propeller. The 
electrical terminals of the winding will be brought outside the motor 
through slip rinss. This winding will hence forth be referred to as 
the shell. 

The rotor ard shell are surrounded by yet another three phase four 


pole winding that contains no magnetic core material. ae. als 





rigidly attached to the motor frame and contains only enough steel to 
provide mechanical support and rigidity. This winding will hence forth 
be referred to as the stator. 

Surrounding the stator is a shield made of a magnetic core material. 
Its purpose is to contain the magnetic fluxes within the motor and 
provide a more efficient return path for the fluxes. If the shield were 
not present the rotating flux field would extend beyond the motor and 
possibly induce hazarious voltages in the surrounding areas which could 
present a hazard both to equipment being operated in the vicinity and 
to mane 

The stator is now surrounded by a water tight enclosure, called the 
pod, that will keep the motor dry and minimize the hydrodynamic drag. 
The magnetic shield could be a part of the pod body. 

In operation, three phase electrical power will be provided to 
the stator winding which will preduce a revolving magnetic field. The 
rotor will lock into synchronism with the stator field and in effect 
become a virtual magnetic core by producing arotating magnetic field 
far greater than could the air core stator by itself. The shell will 
act as a wound rotor induction armature under the influence of the 
rotating magnetic field and drive a super cavitating propeller directly 
as shown in figure (1). Figure (2) shows the geometric configuration 


of the field, damper, shell stator and shield. 
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Figure (2) 
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‘JII MATHEMATICAL MODEL OF THE MOTOR 


The motor is modeled electrically as shown schematically by figure 
(3). Figure (3) shows the axes of the various windings. The rotor has 
a direct axis (in the direction of positive flux lines of the field) ayatil 
a quadrature axis which leads the direct axis by ninety degrees. The 
direct axis is at an angle 96 with respect to the stator ‘A' phase 
axis as shown in figure (3). 

The shell 'A‘' phase axis is at an engle dg with respect to the stator 
‘A’ phase axis. 

The various mutual inductances are assumed to be sinusoids dependent 
upon the angular orientation of the axes involved. The mutual inductances 
are maximum when the axes are aligned, zero when they are ninety degrees 
apart, and negative (minimum) when they are one hundred-eighty degrees 
apart. Equation set (1) gives the flux linkages for the various circuits 
when the currents in all circuits are considered. Equation set (2) gives 
the voltages at the terminals of the various circuits when all currents in 
all circuits are accounted for and the circuit relation Vy = PA, +i R. I, 
is utilized. The symbol p is a operator representing the time derivative 
d/dt. All other symbols utilized are idetified in the Glossary of Terms, 
Mele (1). 

The rather formidable equation sets (1) and (2) can be significantly 
simplified by performing a Park's transformation on theme The Park's 
transformation equations are given here as equation set (3). When the 
transformation of equation set (3) is applied to equation sets (1) and 


(2), equation sets (5) and (6) result. 
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Table I Glossary of Terms 


Note: where a dual assignment of usage is indicated, the context should 


make it clear to the reader the intended usage of the symbols. 


~ynbols 


I Current 
L Self Inductance 
1 length 


M Mutual Inductance 


N Coil Turns/ angle as defined in figure (3) 
iP Permeance 

p Differential Operator d/dt 

R Resistance/ Radius 

S slip 

Jb Torque 

a time 

U Dummy Variable used in the Park's Transformation Equation 
V Voltage 

W Angular Velocity 

Xx Reactance/ Ratio of Reactances 

6 Angle Defined in Figure (3) 

d Flux 

TI tn on 


gd Angle Defined in Figure (3) 


yy Per-Unitized Flux 





Table I continued 
wubscripts 

a Phase a/ Stator Armature 
b Phase b/ Base Quantity 

c Phase c 

d Direct Axis 

£ Field 

2 Arbitrary Index 

5 BeTined an Text 

k Wet nNedma bcs 

a Detinedgim Text 

m Defined in Text 

o) Zero Sequence/ Synchronous Frequency/ Defined in Text 
P Per-Unitized Quantity 

ic Rotor 

Ss Shell 
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Zi 
Equation set (4) can now be obtained directly from equation sets (5) 


and (6). 

A further simplification may be made to equation set (5) by 
per unitizing it in accordance with reference (3). First, equation set (5) 
is re-written as equation sets (7), (8) and (9). 

Then, the ordinary variables in equation sets (7), (8) and (9) are 
normalized in accordance with reference (3) by dividing each variable 
quantity by a corresponding base quantity. The results are shown as 
equation sets (10), (11) and (12). 

It is highly desirable that the resulting per-unitized mutual 
inductances exhibit the property of reciprocity, that is Meg = Mase 
This imposes the following restrictions on the base quantity selection. 


Hf if 


a 5 = = 
E 2 oa Tp Voss ne] 0 View fb 


2 
2 ens ab 


Additionally, the following conditions on the base quantities are 


' arbitrarily imposed to simplify the per-unitized equations. 


a = 3 Mas Wo Fst bo Mar “o trp = Mar Yo TE = S Te bee Wo ab 
am 2 V V V V 
ab ab ab ab 


Thus, at this point, six constraints have been pieced on the eight base 
quantities. This leaves only two base quantity parameters free to be 
chosen arbitrarily. Once these two more conditions are specified the base 
quantities will be completely determined. Further consideration of base 
quantities is given in appendix (H). 

Equation sets (10), (11) and (12) may now be written as equations sets 


(13), (14) and (15). 
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Appendix (A) shows that hop eae and X= X 


So equati 
ee sf quation 


ri 
sets (13), (14) and (15) can be re-written as equation sets (16), (17) 
and (18) respectively. 

If the base quantities utilized to obtain equation sets (10), (11) 
and (12) are applied to equation set (4), equation set (19) results. 

Now, if equation set (19) is solved for the time derivatives of 
the fluxes, equation set (20) results. 

Equation sets (16), (17) and (18) can be inverted so that fluxes 
are the dependent variables and currents are the independent variables. 
The resulting equations are shown as equation sets (21), (22) and (23); 
where, x. ( eae Xe ) 


Sci X aCe Xo) 
2, a S 


Coax. ) 


Equation sets (21), (22), (23) and (20) can be combined to obtain 
equation set (24). 
The direct and quadrature axis voltages of the stator armature are 


given by: ages = Vv, cos (@) and UeeD 2 = VS (cy) 
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The dynamic mechanical equations of the shell and rotor in per 


a — 
unit form are: p( 6) Pa ee ie H,, 
Z = i 
and p ( 0=-N) W (75> T, /2 He 


Appendix (B) provides the following expressions: 


ee) c = 
Tp W (2 * sq eed ee 
Ko 
top ; a | sd = fF ad ese ) a Tp 


These equations along with the equations of equation set (24) form 
a complete set of differential equations that mathematically describe 
the motor's model. It is assumed in this thesis that all loadings and 
transient states will be symetricel and balanced. Thus, W and 


a0 
e will be assumed to be zero. Additionally, & . will be considered 


£ 
a constant. This being the case, the given differential equations. 
reduce to a ninth order non-linear set of differential equations. 
These equations can be solved once all the constant parameters are 
known. This solution is rather straight forward for the steady state 
condition, although a computer minimizes the computational effort. 
The transient state defies a hand solution. The digital computer 
can be used to simulate the transient state once a set of initial 
conditions has been determined. 


References (8), (9), (10) and (12) provided valuable background 


material that aided in the development of the equation in this section. 
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IV Determination of Machine Parameters 
The size of the motor is determined in Appendix (C) and the circuit 
parameters are determined in Appendix (D). The results are summarized in 


Table (2). 





Table 2 Machine Size and Circuit Parameters 


Inner Radius of Field 


Outer Radius of Field 


Average Radius of Damper 


Inner Radius 
Outer Radius 
Inner Radius 
Outer Radius 
Inner Radius 


Outer Radius 


of Shell 
of Shell 
of Stator 
of Stator 
of Shield 


of Shield 


Electrical Length 


Rap 


R 
Sp 
R 
rp 
“oR 
A. 
J 


4.0 inches 
6.5 inches 
7-0 inches 
7°5 inches 
10.5 inches 
11.0 inches 
13.5 inches 
14.5 inches 
20.0 inches 
74 feet 
0.20 
0. 588 
26.87 
0.0 (super-conducting winding) 
70799 
3.685 
4.685 
10.08 
0.738 seconds 


0.0738 seconds 





- 
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V Steady State Characteristics 


When the motor is running in steady state, the governing 
differential equations of equation set ( 24) can be reduced by setting 
all time derivatives of the fluxes to zero. Additionally, the torque 
on the rotor must be zero and this directly implies era = W on i QO. 


q 
Slip, s, is defined defined in the conventional manner: that is, 


ng S 
: lem where W = p(@)= We and Wo = p(¢) 
wo - PCS) 
in the steady state, then s = = 7 ne As can be seen from 
O 


figure (3), 02=N+¢ so PCRs) = We -~p(¢), so finally p(N)=s Woe 
Taking equation set ( 24) and setting the time derivatives of 


fluxes to zero, SR ie = P va = 0, and using the relations 


pN = s Wo » pod = Wo and I p = Yep i) Rep » equation set (25) results. 


f 
The torque equation of the shell reduces to, 
Xo 
Te ee le : 
sp oe sd aq and by using the fourth equation of 


equation set (25), the expression for torque can be written as, 


R X Ww 
Ss W 


sp “3 
Wa /s qT. 

Solving for yields oad = —— “tba. 
aq aq He Secs 


Assuming that slip, s, stator torque, Tey? and field current, Lis 
are the independent variables, then equation set ( 25) can be solved 
for the fluxes and voltages which result in equation set ( 26). 


The digital computer is a valuable tool in evaluating the expressions 


of equations set ( 26) numerically. Equation set (26) and the machine 
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parameters listed in Table (2) along with the load characteristics 
determined in appendix (E) were injected into an IBM model 370 
computer use the WATFIV program listed in appendix (CGC). 


The results of the computer analysis are shown in the results 


section as figure (4) through (15). 





VI Transient State Characteristics - 

Equation set ( 24), table (2) and the shell and rotor torque 
expressions of Appendix ( B) provide the mathematical discription 
of the transient state. Given these equations and a complete set of 
initial conditions and driving functions, the system is deterministic. 
The initial conditions are determined in Appendix (F). 

The transient state was simulated in four distinct evolutions: 

Meeenaeor start up with ‘ee = infinity 

Ze Shell start up with rotor at full speed and at various values of 

Begs ‘ 

Sic Crash back. 

4, Three phase fault at the stator armature terminals. 

The simulations were all conducted using an IBM model 370 digital 
computer and the IBM supplied program CSMP (Continuous Systems Modeling 
Program). The program listings are shown in Appendix (G). | 

The rotor start up simulation assumed that the shell circuit is 
open and in effect not there. Thistas equivalent to starting a 
conventional single armature machine. The results of this start up 
ar shown in figure (16). 

The shell start up simulation was conducted by first having the 
rotor running in steady state then placing various values of finite 
resistance in the shell circuit. The results of this simulation are 
shown in figure (17). 

The crash back simulation first assumed a steady state full power 


condition for the motor, then the phase sequence of the stator terminal 
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voltage were reversed. Tne results of the crash back simulation are 
shown in figure (18). 

The three phase fault Simulation again assumes a steady state full 
power condition for the notor, then the three terminal leads of the stator 
armature are electrically shorted. In effect ie =V.,.=V 0. 


ab ac 


The results of the three phase fault are shown in figure (19). 





VII Results 
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Figure (4) through figure (15) show the steady state analysis 


results. In these figures, s, I 


Figures (4), 
Figures (7), 
Figures (10), 


Figures (13), 


and Tare the independent variables. 
fp Sp 

(5) and (6) show the terminal voltages. 

(8) and (9) show the terminal currents. 


(141) and (12) show the power factor. 


(14) and (15) show the rotor angle. 


Figure (16) through figure (19) show the angular velocities of the 


shell in the transient 
Figure (16) 
Figure (17) 
Figure (18 ) 


Figure (19 ) 


State. 

shows the rotor start up transient. 
Shows the stator start up transient. 
shows the crash back transient. 


shows the three phase fault transient. 
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VIII DISCUSSION OF RESULTS AND CONCLUSIONS 


Figures (4) throuch (15) demonstrate that the dual armature 
motor can be characterized by a set of steady state condition curves. 
These curves can assist in the rational selection of an operating 
condition for the motor. The data supplied by these curves can be 
cut many different ways; for instance, Lop ORE Tp could be considered 


as dependent variables and un or i, could be treated as independent 
variables. In cutting the data different ways additional data points 
may be desirable, in this case the WATFIV program in Appendix (G) could 
be used if the appropriate parameters were adjusted. 

Figures (16) through (19) provide the transient state results. 
Figure (16) shows that the rotor on start up under goes large angular 
velocity excursions with the attendant large acceleration forces. At 
One point the rotor actually runs tackvardse Eventually the rotor 
settles down to its mechanical synchronous speed of 1800 RPM. 

Figure (17) shows that for the shell start up there iS some critical 
parameter Ss between 10 pu and 20 pu where the shell starts up ina 
Stable fashion at the fastest possible rate. For values of Boas less than 
this critical value, the rotor tumbles or pulls out of synchronism with the 
Stator field and results in a loss of starting torque and the introduction 
of a pulsating torque. For Re» = 0.588 pu (its minimum value) the motor 
rotor does not start at all but oscillates about the zero point. 

Figure (18) shows the crash back Wop curve. The results of this 
curve beyond 1.8 seconds are questionable. The computer simulation appears 


to have gone astray at this point, perhaps due to too large of a step 


meze on the Simulation integration routine or due to the proximity of W_ 





to zero. The rotor and shell both move at approximately the same velocity 
in this maneuver. 

Figure (19) shows the three phase fault at the stator armature 
terminals. Again the computer simulation appears to go astray at about 
3.7 seconds so results beyond that point are questionable. 

The preliminary results of this thesis indicate that the dual 
armature super conducting motor is indeed a candidate for hydrofoil 


propulsion and warrants further investigation. 





IX RECOMMENDATIONS 


Additional areas that could be investigated using the equations 
developed in this thesis may be: 

i, Determine steady state conditions as a function of Ne » The 
present investigation considered Ss = 0,588 pu in the steady state 
condition. 

Ze Write a computer program that will accept machine dimensions as 
input and return machine parameters for use in the steady state 

or transient analysis. 

36 Instead of holding © ep constant as done in this thesis in the 
transient state investigation, hold Lep constant. 

L. Investigate current densities in the various windings during 
the transient periods to see if they exceed tolerable levels. Devise 
a scheme for controlling excess currents during transients. 

5e Determine the static stability limit of rotor angle. 

6. Refine load characteristics. 

oP Vary machine dimensions to optimize the design. Criteria for 
optimizing could also be developed as a side issue. 

8. Introdue hydrofoil dynamics into the model. 

96 Investigate speed control by varying No the electrical fre- 
@uency Of the; prime mover. 

Most of these subject areas would apply at an under graduate project 


level. 
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Appendix A Shell Constraints 


Equation sets (13), (14) and (15) contain parameters that are 
not independente The following development assumes that flux 
leakages between adjacent circuits is negligible. 

Assuming that the stator, shell and field axes are aligned so as 
to eliminate the sinusoidal dependency of the mutual inductances, the 
following relations can be written. 


sa eon - 
( L, Mo ) 2 Ns C Poa 


a ee 
a we ME ( Feacer ) 


ayy SCE? CRRA) 
Meee No et ODF) 
oe AE, BC a AS) 
hag NE. N. ( Eec.F a) 
ages Se N. ( rot) 


Nica Nea TO ag 


The 2 factor appearing in the first equation is characteristic of 
three phase windings. 

Now, if the various leakage permeances may be considered negligible 
and are set to zero, the above equations may be used to derive the 


following relations. 
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These equations may now be per=-unitized as follows: 
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Appendix B Derivation of Torque Expressions 


The torque on any winding in per unit is given by 


ee, {445 aed tigp 


Applying this to the stator and shell windings and substituting the 
expressions for current from equation Sets (21) and (22) and canceling 


terms where possible yields: 





ie an 
Tap ie (2 34 Eq ~ eer ale a 
ao 
‘sp it, Ceca aon 7 “are Pa ) e395 
meso equilibrium of forces requires T ¢+T +T =O, 
ap SD ine 


Bolving the above equations for the torque on the rotor gives 


X 
O 
Tp Hs © sq a as 7 © a = vg )s 








Appendix C Motor Sizing is 

References (11) and (14) contain equations for the power rating of 
an air core, three phase synchronous motor, They can be simplified and 
adapted to the present application by assuming the followings 

Ae number of pole pairs = 2 

be winding angle of field = 120° 

Ce Winding angle of shell = 60° 

d. winding angle of stator = 60° 

Ce all winding electrical lengths= 1 


ts the following parameters are defined as 


x = Bea Roe , yr Rey / Bee , and z= R,,/R, » The radius 
dimensions are shown in figure (C-i). 
The resulting power rating expression is, 
P= 12 65407 J*_ i. C.'R (18 hy it 
V Om ab i 22 s fo E 


Ly R 
where C, = = ( -~Inx+t (ne) | ~22 )) , 
sh 


2 


Ve / 2 
cae ake SIN(Y) =./1 - Xx, cos (Yi); 
COS(Y) = power factor , 
x, =, /0.375 etatom “ee ( 320 9! 
= e == == o 
a Jp C., 1- Rag 


by 4.2 R 
mae Cees ints ee co ty, 
S Ly 8 Ron 
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Figure (C-1 ) 
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The current densities were chosen based on the values used in 
reference (11) and are; J. 2 2.50 x 10° amps /meter* 4 


: amps /meter~ 3 


and J, = 1.25 x 10 
Synchronous electric frequency was chosen to be Wo = 377 radians/sec.e 
Power factor was chosento be PF = 0.8 .o 
The power was chosen based on the requirements specified in reference 


(4) and is P= 21,000 horespower per motor. 


The radii selected are; 


Res = 4,0 inches Reo = 6.5 inches 
Ros = 7.5 inches Reo = 10.5 inches 

R_. = 11.0 inches R__ = 13.5 inches 
ai 3.0 

Ry = 7.0 inches Ran = 14.5 inches 


Finally, if the above assumptions are made and the power rating 


expression is solved for 1, the result is 1 = 7.4 feet. 
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_ Appendix D Circuit and Mechanical Parameters of Motor 


References (11) and (14) contain equations for determining the 
inductances (both self and mutual) for an air core, three phase, 
Synchronous motor. These equations are infinite series Summations 3 
however, the lead term, n=1, dominates and all other terms are small, 
These equations can be simplified and adapted to the present application 
by assuming, in addition to the assumptions of appendix (C), that: 

ae the leading term of the infinite series dominates, 

De the shield is made of magnetic iron and 

Ce the theory of superposition holds, 

The resulting inductance expressions are listedinTable (D -1), 
The dimensions are the same as discussed in appendix (C) and shown in 
figure (C~-1). If these dimensions are now used to solve the expressions 
in table (D-1), table (D-2) results. 

If the expressions of table (D=2) are per unitized using the 
base quantities of appendix (H), table (D-3) results. 

The definitions of oe Xs Xo» x and KX, given on page 25 can now be 
evaluated and are; Ks = 7.799 

xX, = 3.685 
Xx, = 14.685 
Xo = 10.08 
Reference (14) provides an equation for R, and is 


Tiel 2 7 
1 


. mp 


————— 
r 4 0. Ry 


where Oo. = 90 t 


q? ty being the damper shell thickness and 0 being the 
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Table D-t Inductances (cont.) 
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Table D=3 


Inductances 


= } old 2 
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conductivity of copper. Assuming GO = 3.6 x 10° mhos/ eter and 
t, = Oel1 inches, R, = 4.42 x 107? nN 1. Per unitizine usinel the 
base quantities of appendix (H) gives, Ray = 26.87 pu. 
In general, the resistance of a wound three phase winding is 
; 2 
given by dy N. 1, 
®eld 
oi shedls wei si i 
assuming the conductivity is 0 = 6x 10° whos /meter and the space factor 
Ke = 0.27 for both the shell and the stator armature, the per unitized 
resistances ares R = 0.200 pu and R  =0.58 . 
ap Sp 
The resistance of the field winding is of course zero because 
it is made of a superconductor, 
An initial estimate of the moments of inertia for the rotor and 
for the shell with shaft, propeller and entrained water included were; 
Ag = 120 slug et" and oe = 12 slug ft a 


For H defined in the conventional manner, that is, 


Zz 
i J. We 
i 2 P. 
H. = 0.738 sec and H = 0.0738 sec. 


Note: the computer listings of appendix (Ge) show Hto be different 
from the above values. This is because Wo was included so the value in 


the listings is 377 times smaller than the actual assumed value. 
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Appendix E Load Characteristics 


The dynamics of the motor and propeller system are 
considered to be uncoupled from the hydrofoil dynamics for the 
time span of the simulations since the time constants of the two 
systems are quite different. Thus the load torque characteristics 
are a function of W. alone e 

For the shell start up,the propeller characteristics are 


assumed to be = z 
i (w/a) ; 


For the crash astern and the three phase fault, the load 


characteristics are assumed as shown below. 


fe rp / 
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Appendix F Determination of Initial Conditions 


An operating point was selected for the full power steady state 
from figures (4) through (15) as s=0.0015, Tp 71-0 pu, and Te pate0 pu. 
Then using the equations of equation set (26) a complete set of initial 
condition were determined and are listed in the computer listings in 
appendix (G). 

FOreroverss tart Up, 1% was assumed tnat Ta? W = 0.05 W 0-05 
and §=0.0 . The remaining initial conditions were calculated using 
Equation Set (26) and are listed in the computer listings in appendix 
ac). 

Part of the rotor start up computer simulation routine was to 
print out the final values of all the time variables at the end of 
the simulation period. These final values were then utilized as the 


initial conditions for the shell start up. The values are listed in 


the computer listings in appendix (G). 





Appendix G Computer Program listings 
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Appendix H Base Quantity Selections 
As mentioned on page ( 21), two remaining base parameters are free 
to be selectede The bases selected ares 
1) P= 21,000 horsepower and 
2) Lop = 1.0 pu when Jp = ee 10° amps /meter~. 
At this point eight conditions have been specified on the eight 
base quantities so the base system is fully constrained and all base 


quantities can be calculatede They are 3 


= a 5 
VES = 19.36 N, volts ie = 26070 0 amps /N, 
Vv. = 41.65 N_ volts Ti. = 4.479 x 10° am /N 
Sb S Sb ‘ PST Ss 
ae = 5-640 Ne volts l= 1.388 x 10° amps /N, 


os HOBBS : 
Ve, = N, volts Top = 20247 x 10° anps/N, 
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